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A  REVIEW  OF'  ONE  DIMENSIONAL  SHAPED 
CHARGE  THEORY 
PART  1  -  JET  FORMATION 


1 .  INTRODUCTION 


The  shaped  charge  warhead,  based  on  the  Munroe  effect,  has  been  used 
for  the  defeat  of  heavy  armour  since  the  second  World  War,  A  particularly 
comprehensive  account  of  the  historical  development  of  these  weapons  has 
recently  been  given  by  Backofen  [1 ] . 


Shaped  charge  warheads  are  currently  in  use  by  the  Australian  Army 
in  the  66  mm  M72L1A2  ( IAW)  and  84  mm  Carl  Gustaf  anti-tank  weapons,  and  it  is 
anticipated  that  an  advanced  shaped  charge  device,  the  103  mm  MILAN,  which 
incorporates  wave  shaping  for  increased  performance,  will  soon  be  added  to  the 
list.  Also,  it  seems  highly  likely  that  future  torpedo  designs  will 
incorporate  shaped  charge  warheads  for  defeat  of  submarine  pressure  hulls,  as 
in  the  Stingray  torpedo  used  by  the  Royal  Navy, 


MRL  has  been  engaged  on  various  investigations  concerning  shaped 
charge  research,  development  and  evaluation  for  some  years.  Examples  of  this 
work  are  the  design  and  fabrication  of  a  38  mm  copper  lined  shaped  charge  [23. 
which  has  since  become  known  as  the  MRL  Standard  Shaped  Charge,  and  the  recent 
use  of  this  charge  for  a  very  illuminating  investigation  on  the  jet  initiation 
and  disruption  of  bare  and  covered  Comp  B  and  other  explosives  by  M.C.  Chick 
and  D.J.  Hatt  [3], 


Interest  in  the  modelling  of  shaped  charge  performance  has  also  been 
stimulated  by  the  use  at  MRL  during  the  last  few  years  of  several  large  two 
dimensional  (2D)  continuum  mechanics  finite  difference  codes  [45-47],  Of 
these,  one  of  the  most  useful  for  shaped  charge  work  is  the  HELP  code.  This 
has  recently  been  used  by  D.L.  Smith  to  model  the  early  staqes  of  jet 
formation  in  the  MRL  38  mm  shaped  charge,  and  a  preliminary  comparison  with 
flash  radiographs  of  the  collapsing  liner  is  extremely  good  [4] .  Fig.  1  shows 
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a  schematic  illustration  of  the  progressive  stages  of  liner  collapse  and  the 
formation  of  jet  and  slug,  and  Fig.  2  shows  one  of  the  recent  MRL  flash 
radiographs  of  the  collapse  of  the  MRL  38mm  Standard  Shaped  Charge  18.5  ys 
after  initiation  of  detonation. 


Research  on  the  basic  understanding  and  refinement  of  shaped  charge 
devices  is  currently  being  pursued  in  many  weapons  laboratories  around  the 
world,  and  much  interesting  work  has  been  published  in  the  open  literature  and 
in  unclassified  reports  during  the  past  decade.  The  purpose  of  this  review 
is  not  only  to  acquaint  the  reader  with  these  recent  developments,  but  also  to 
provide  a  reasonably  detailed  account  of  the  basic  hydrodynamic  theory  of 
shaped  charges  which  was  published  in  the  late  1940's  and  early  1950's.  It  is 
hoped  that  this  review  will  be  of  interest  to  people  either  using,  designing 
or  modelling  shaped  charges,  as  well  as  providing  sufficient  background 
material  for  the  understanding  and  formulation  of  a  one  dimensional  (ID)  code 
for  the  analysis  of  shaped  charge  behaviour. 


Modelling  of  classical  shaped  charges  by  the  use  of  so  called  ID 
codes  has  become  quite  common  in  the  past  few  years.  Possibly  the  best  known 
example  of  such  an  approach  is  the  BASC  code  used  by  the  Ballistic  Research 
Laboratory  in  the  U.S.A.,  which  has  been  extensively  documented  by  Harrison 
[51.  Other  codes  include  JETFORM,  used  by  RARDE  in  the  U.K.  [6],  DESC-1 , 
designed  by  Car leone  et  al  of  Dyna  East  Corp  [7] ,  and  the  TB/ISL  code  of 
Hennequin  in  France  [8].  As  pointed  out  by  Carleone  et  al  [9],  these  codes 
are  an  invaluable  addition  to  the  overall  strategy  for  the  design  of  new 
shaped  charge  devices.  Whilst  the  2D  codes  are  capable  of  providing  very 
detailed  information,  they  suffer  from  the  disadvantage  of  reguiring  large 
computers  and  long  run  times,  and  as  a  consequence  are  not  particularly  well 
suited  to  extensive  parametric  studies.  By  contrast,  the  ID  codes  evaluate 
simplified  analytical  expressions  rather  than  a  fully  fledged  finite 
difference  scheme,  which  the  name  might  suggest.  ID  codes  can  be  stored  on 
the  smallest  of  computers,  require  only  seconds  run  time,  and  are  ideal  for 
parametric  studies.  An  application  of  the  ideas  discussed  in  this  report  to 
the  MRL  Standard  Shaped  Charge  is  currently  in  progress,  and  will  be  reported 
shortly  [10]. 


An  outline  of  this  report  is  as  follows:  section  2  describes  the 
basic  steady  state  hydrodynamic  theory  of  jet  formation  published  by  Birkhoff, 
MacDougall,  Pugh  and  Taylor  in  1948  [11).  The  description  here  is  rather 
detailed  as  this  paper  provides  the  basic  understanding  of  shaped  charge 
behaviour.  The  theory  was  modified  in  1952  by  Pugh,  Eichelberger  and  Rostoker 
to  describe  the  stretching  of  the  jet  in  free  flight,  and  this  resulted  in  the 
"non-steady  state"  theory  of  jet  formation,  commonly  referred  to  as  the  PER 
theory  [12].  This  work  is  described  in  section  3.  Before  the  PER  theory  can 
be  utilised,  expressions  must  be  found  for  the  liner  velocity  once  it  has  been 
set  in  motion  by  the  detonation,  and  a  variety  of  methods  for  achieving  this 
are  described  in  section  4.  The  velocity  of  the  jet  tip  is  an  important 
parameter  in  shaped  charge  work,  and  its  calculation  is  complicated  by  an 
inverse  velocity  gradient  within  the  leading  portion  of  the  jet,  which  causes 
the  compression  of  the  jet  and  the  formation  of  a  compact  jet  tip  particle. 

The  treatment  of  this  effect  is  described  in  section  5.  Finally,  the  early 
theories  of  Birkoff  et  al  and  Pugh  et  al  assumed  that  the  flow  was 
incompressible,  in  which  case  a  coherent  jet  must  always  be  formed.  The 


modification  to  this  result  obtained  by  allowing  the  flow  to  be  compressible 
is  described  in  section  6. 


This  review  is  devoted  entirely  to  the  problem  of  jet  formation,  and 
no  attempt  has  been  made  to  include  recent  work  on  the  important  topics  of  jet 
breakup  and  jet  penetration.  It  is  intended  that  a  review  of  such  material 
will  appear  in  a  future  MRL  report. 


2.  STEADY  STATE  THEORY 


As  mentioned  in  the  introduction,  the  first  fairly  complete 
mathematical  theory  for  a  classical  shaped  charge  was  published  in  the  open 
literature  in  1948  by  Birkhoff,  MacDougall,  Pugh  and  Taylor  [11].  Their 
analysis  was  based  on  the  following  assumptions: 

(i)  The  pressure  on  the  liner  due  to  the  explosive  is  so  great  that  the 
yield  strength  of  the  metal  is  negligible,  hence  the  metal  will 
behave  hydrodynamically.  (Note  that  this  does  not  imply  that  the 
metal  is  a  liquid  during  the  formation  process). 

(ii)  The  metal  flow  is  non-viscous  and  incompressible. 

(iii)  The  pressure  on  the  collapsing  liner  does  not  increase  the  length  of 
any  portion  of  the  liner. 

v iv)  The  impulse  from  the  detonation  wave  acts  instantaneously,  so  that 

the  acceleration  of  the  liner  can  be  neglected. 

(v)  Each  element  of  the  liner  collapses  with  a  constant  velocity^ 

which  is  independent  of  its  original  position  along  the  cone  axis. 

Several  of  these  assumptions  have  been  modified  or  relaxed  in  more 
recent  refinements  of  the  theory,  and  these  will  be  discussed  in  detail  in 
later  sections  of  this  report.  For  the  moment,  we  consider  the  original 
theory  as  presented  by  Birkhoff  et  al. 


Consider  Fig.  3  (adapted  from  reference  [3]),  which  shows  the 
detonation  wave,  with  velocity  traversing  the  length  of  the  cone.  PN 
represents  the  velocity  of  the  liner.  From  assumption  (v)  we  know  that 
each  element  of  the  liner  between  A  and  P  (the  current  position  of  the 
detonation  front)  has  been  compelled  to  move  with  the  constant  velocity  _VQ. 
(ie_PN[=  AA ' ),  and  from  assumption  (iii)  we  know  that  the  length  AP  is  eoual 
to  the  length  A’ P.  We  now  show  that  PN  bisects  the  angle  between  the 
uncollapsed  liner  PB  and  the  collapsing  liner  PM.  As  PN  is  parallel  to  AA' 
we  have  angle  BPN  =  angle  PAA' ,  and  angle  NPA'  =  angle  PA' A.  Also,  as  AP  = 
A'P,  then  angle  PAA'  =  angle  PA' A,  ie  angle  BPN  =  angle  NPA',  which  was  to  be 


Having  found  the  direction  of  we  now  derive  an  expression  for 
its  magnitude,  again  using  assumptions  (iii)  and  (v).  Consider  Fig.  4,  which 
shows  the  detonation  wave  having  advanced  from  p  to  p'  in  the  time  ta/un  for 
the  liner  segment  launched  from  p  to  hit  the  axis  at  N.  Construct  the  line 
P'N.  Again  by  assumption  (iii),  we  have  P'P  =  p'N.  Also  construct  the  line 
PL,  which  is  perpendicular  to  the  original  liner  surface  at  P,  and  the  line 
P'K,  the  perpendicular  bisector  of  PN.  Define  angle  NPL  =  6^  and  PP'K  =  ^2* 

It  is  easy  to  show,  by  summing  angles  at  P'  and  p,  that  ^  =  $2  =  V2  (@-°0/ 
where  0  is  the  angle  which  the  collapsing  liner  makes  with  the  axis,  and  a  is 
the  original  cone  half  angle.  From  triangle  P'PK  we  then  have 

sin  6  =  V  /(2U) ,  (2.1  ) 

o 


where 


U  =  U  /cos  a  . 
D 


Equation  (2.1)  is  fundamental  to  shaped  charge  work  and  was  first  derived  by 
G.I.  Taylor  [14],  and  6  has  since  become  known  as  the  Taylor  angle.  It  has 
recently  been  extended  by  Randers-Pehrson  [15],  and  by  Chou  et  al  [16],  to 
apply  to  non  steady  state  conditions.  These  corrections  will  be  discussed  in 
detail  in  the  next  section. 


Next  we  derive  expressions  for  the  velocities  of  the  jet  and  slug. 

We  denote  the  velocity  of  the  collision  point  M  by  ,  and  then  move  to  a 
frame  of  reference  in  which  M  is  at  rest  by  imposing  a  uniform  translational 
velocity  of  -  to  m  this  frame  each  segment  of  the  liner  PM  then  has 

a  velocity^  and  moves  along  the  line  PM.  Fig.  5  illustrates  the  velocity 
relationship.  Applying  the  sine  rule  we  have 

V  V„  V, 

o _ 2 _  _  _ 1 _ 

sin  ®  sin  [-^  -  (c  +  6)]  sin  [~  +  (a  +  6  -  0)  ] 

from  which  we  find 

V.  =  V  cos  (a  +  5)/sin  0  (2.2) 

2  o 

V  =  V  cos  (0  -  a  -  5) /sin  0  (2.3) 

1  o 

Because  of  assumption  (v),  an  observer  in  the  frame  of  reference 
moving  with  the  velocity  will  see  the  steady  state  situation  depicted  in 
Fig.  6,  ie.  two  fluid  jets  with  velocities  colliding  at  an  angle  20.  This 
situation  has  been  described  by  Milne-Thomson  [17],  and  the  result  of  the 
collision  will  be  two  fluid  streams,  a  "jet"  moving  to  the  right  with  velocity 
^2,  and  a  "slug"  moving  to  the  left,  also  with  velocity  J/,,.  Hence,  in  the 
stationary  system  of  coordinates,  the  jet  will  be  moving  to  the  right  with  a 
velocity  given  by 


Vj  " 


V,  +  V, 


(2.4) 


while  the  slug  will  also  be  moving  to  the  right  with  the  slower  velocity 


(2.5) 


V  3  V  •  V 
vs  V1  v2 

Using  equations  (2.2)  and  (2.3)  in  (2.4)  and  (2.5)  we  find 

V.  =  V  cos  (a  +  5  -  0/2)/sin  (0/2),  (2.6) 

3  o 

V  «  V  oJn  (a  +  6  0/2)/cos  (0/2).  (2.7) 

s  o 

For  the  simple  case  treated  here  in  which  is  independent  of 
position  along  the  liner  we  have  already  found 

6  =  V2  (0  -  a),  (2.8) 

in  wiij.ch  case  equations  (2.6)  and  ' 2.7)  take  the  simple  form 


V..  =  Vq  cc-i  (a/2)/sin  (0/2),  (2.9) 

V  =  V  sin  (a/2)/cos  (0/2).  (2.10) 

so 


The  division  of  liner  mass  between  jet  and  slug  is  governed  by  the 
conservation  of  momentum.  Let  m  be  the  liner  mass  per  unit  length 
approachira  the  junction.  Let  nu  be  that  part  of  m  going  into  the  jet  and  mg 
be  that  goiij  into  the  slug.  By  equating  the  horizontal  components  of 
momentum  before  and  after  passing  the  stagnation  point  M  in  the  moving  frame 
we  have 


(2.11 ) 

(2.1  2) 


nij/m  =  sin2{0/2),  (2.13) 

m  /m  =  cos2(0/2).  (2.14) 

s 

It  is  important  to  point  out  that  this  steady  state  analysis  is 
strictly  true  only  for  the  case  of  a  wedge  shaped  liner.  Consider  the  cross- 
section  of  the  charge  shown  in  either  Fig.  3  or  Fig.  4.  The  wedge  shaped,  or 
linear  charge,  is  defined  by  supposing  the  cross-section  to  be  the  same  for  an 
infinite  distance  perpendicular  to  the  cross-section,  while  the  conical  case 
is  obtained  by  rotating  the  cross-section  around  the  line  of  symmetry.  in 
order  for  the  whole  of  the  collapse  process  to  appear  stationary  and  the 
simple  situation  depicted  in  Fig.  6  to  be  applicable  it  is  necessary  for  the 
total  mass  per  unit  distance  along  the  axis  to  be  constant.  This  is  true  for 
a  wedge  shaped  charge,  but  is  only  approximately  true  for  a  conical  liner  of 
constant  thickness.  Nevertheless,  the  application  of  this  theory  to  conical 
liners  is  remarkably  accurate,  as  we  will  see  in  later  sections. 
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m  summary,  the  steady  state  model  predicts  that  both  the  velocities 
and  cross-sectional  thickness  of  the  jet  and  slug  will  be  constant. 
Experimentally  however  this  is  known  not  to  be  the  case:  the  jet  in  particular 
has  a  marked  velocity  gradient  along  its  length,  the  tip  of  the  jet  travelling 
some  four  to  five  times  faster  than  the  base,  which  leads  to  considerable 
elongation  of  the  jet  in  flight  before  target  penetration.  This  fact  was  known 
by  Birkhoff  et  al  and  was  used  by  them  in  their  description  of  jet 
penetration,  but  was  not  included  in  their  theory  of  jet  formation.  The 
modification  to  the  simple  steady  state  theory  to  account  for  the  gradient  in 
jet  velocity  was  made  by  Pugh,  Eichelberger  and  Rostoker  [12]  and  is  described 
in  the  next  section. 


3.  NON  STEADY  STATE  THEORY 


The  velocity  gradient  within  the  jet  was  satisfactorily  explained  by 
Pugh  et  al  [12]  by  modifying  assumption  (v)  as  follows: 


The  velocities  with  which  the  various  elements  of  the  cone  liner 
collapse  when  they  are  struck  by  the  detonation  wave  depends  upon  the  original 
position  of  the  element  in  the  cone. 


This  is  certainly  a  reasonable  assumption  to  make  as  the  ratio  of 
the  explosive  mass  to  liner  mass  decreases  continuously  from  apex  to  base  and, 
bearing  the  Gurney  formulae  in  mind  [18],  we  would  expect  a  corresponding 
decrease  in  liner  velocity  (remember  that  we  are  still  assuming  that  the 
period  of  acceleration  is  negligible) .  The  decrease  in  velocity  of  each 
segment  hitting  the  axis  results  in  a  corresponding  decrease  in  velocity  of 
each  jet  element  as  it  is  formed,  hence  the  velocity  gradient  within  the  jet 
is  explained. 


Because  VQ  is  a  function  of  position  along  the  cone  axis  it  is  now 
impossible  to  find  a  set  of  coordinates  in  which  the  collapse  appears  as  a 
true  steady  state  process.  This  means,  in  particular,  that  it  is  no  longer 
possible  to  use  Bernoulli's  theorem  to  simplify  the  hydrodynamics,  as 
discussed  by  Birkhoff  et  al,  and  which  is  fundamental  to  the  simple 
description  of  the  collapse  as  pictured  in  Fig.  6.  Pugh  et  al  overcome  this 
problem  by  using  the  same  approach  as  Birkhoff  et  al,  but  applying  it  to  a 
succession  of  individual  zonal  elements.  Each  element  is  then  considered 
independently  in  the  appropriate  constant  velocity  coordinate  system.  As  a 
result,  equations  (2.6)  and  (2.7)  for  the  velocities  of  jet  and  slug  remain 
valid  in  the  non  steady  state  theory,  where  now  VQ  is  a  monotonic  decreasing 
function  of  x,  where  x  measures  distance  along  the  cone  axis.  Equations 
(2.13)  and  (2.14)  for  the  distribution  of  mass  between  jet  and  slug  also 
remain  valid  in  the  new  theory. 


Another  consequence  of  the  decrease  in  liner  velocity  with  distance 
along  the  axis  is  that  the  collapse  angle  6  is  no  longer  constant,  but 
increases  as  the  collapse  proceeds.  This  is  easily  explained  by  referring 
again  to  Fig.  4.  Consider  the  element  at  Q,  if  this  collapsed  with  the  same 


velocity  as  the  element  from  P,  then  it  would  reach  S  when  the  element  from  P 
reached  N,  and  the  collapsing  segment  of  liner  NSP'  would  form  a  straight 
line.  However,  as  the  collapse  velocity  at  Q  is  less  than  that  at  P  the 
collapsing  liner  adopts  the  shape  NRP' ,  which  means  that  the  angle  0  which  the 
liner  makes  with  the  axis  at  N  is  greater  than  the  angle  0  which  it  would 
have  made  with  the  axis  if  the  collapse  velocity  were  constant. 


To  calculate  the  variation  of  0  with  distance  along  the  cone 
Pugh  et  al  adopt  the  following  procedure.  Let  the  cylindrical  coordinates  of 
R  in  Fig.  4  be  (r,  z)  and  the  coordinates  of  the  original  position  Q  of  R  in 
the  liner  be  (x  tan  a,  x) .  The  connection  between  the  two  sets  of 
coordinates  is  given  by 


z  =  x  +  V  (t  -  T)  sin  A  (3.1) 

r  =  x  tan  a  -  VQ  (t  -  T)  cos  A  (3.2) 


where  t  is  the  time  (t  =  o  corresponds  to  the  detonation  wave  passing  the  apex 
of  the  cone),  T  =  x/U  ,  and  A  *  a  +  5  (x).  The  slope  of  the  collapsing  liner 
is  given  by  3r/3z,  and  the  collapse  angle  0  is  defined  by 


tan  0 


(If) 


izJ  r=o 


(3.3) 


From  equations  (3.1)  and  (3.2)  we  have 


.  tan  a  -  V  ' ( t  -  T)  cos  A  +  V  cos  A/U  +  V  A'  ( t  -  T)  sin  A 

3  r  _  _ o _ o _ D  o _ 

3z  “  1  +  V  '  (t  -  T)  sin  A  -  V  sin  A/0  +  V  A'  (t  -  T)  cos  A 

o  o  D  o 


where  the  prime  denotes  differentiation  with  respect  to  x.  The  condition 
r  =  o  is  equivalent  to 


t  -  T  =  x  tan  a/V  cos  A  .  (3.5) 


Substituting  (3.5)  and  (3.4)  into  (3,3)  we  find 

sin  0+  -  x  sin  a  (1  -  tan  A  tan  5)  V  '  /V 

tan  0  =  - - - - - -  f  3 .6) 

cos  0  +  x  sin  a  (tan  A  +  tan  6)  V  '/V 

o  o 

-f 

where  0  =  a  +  26  and  is  equal  to  the  value  of  the  collapse  angle  in  the 

steady  state  case.  From  the  discussion  already  given  we  know  that  vQ  will  be 
negative,  while  the  expression  in  the  brackets  in  equation  (3.6)  will  be 
positive  for  usual  shaped  charge  dimensions,  hence  equation  (3.6)  shows 
that  0  will  increase  as  x  increases. 


It  should  be  noted  here  that  the  analysis  presented  above  applies 
only  to  the  case  of  a  plane  detonation  wave.  Generalizations  of  these 
results  have  been  made  by  Allison  and  Vi tali  [19]  for  point  detonation,  by  the 
present  author  for  a  spherically  converging  detonation  wave  [20],  ..nd  by 
Behrmann  for  an  arbitrarily  shaped  detonation  wave  and  liner  [21]. 

An  interesting  development  of  equation  (3.6)  has  also  recently  been 
made  by  Hirsch  [22] .  Using  simple  trigonometric  identities  he  has  shown  that 
the  equation  for  tan  6  can  be  written  in  the  form 

tan  6  =  tan  (3+  +  A3),  (3.7) 

where  A3  is  defined  by 


tan  A3  =  f-  x  sin  a/ (cos  A  cos  6))(V  '/V  ). 

>  '  O  O 


(3.8) 


This  simple  development  clearly  shows  how  the  velocity  gradient  of  the  liner 
continuously  increases  the  collapse  angle  3  to  values  greater  than  the  steady 
state  value  3+.  The  introduction  of  the  angle  A3  also  simplifies  the 
expression  for  the  velocity  of  a  liner  segment  in  the  coordinate  system  moving 
with  the  collapse  point  (equation  2.3),  and  when  this  is  considered  with  the 
criterion  for  coherent  jet  formation  (which  will  be  discussed  in  section  6)  it 
allows  Hirsch  to  make  some  interesting  observations  on  the  effect  of 
confinement  on  the  jet  tip  origin  and  its  variation  with  o. 


Hie  derivation  of  the  Taylor  angle  formula,  eq  (2.1),  clearly  shows 
that  this  relation  is  valid  only  under  steady-state  conditions.  An  extension 
of  this  formula  to  non  steady  conditions  has  been  made  empirically  by  Randers- 
Pehrson  [15],  and  analytically  by  Chou  et  al  [16].  Both  approaches  take  into 
account  not  only  the  gradient  in  liner  velocity  along  the  length  of  the  liner, 
but  also  the  finite  acceleration  of  each  liner  segment. 


Randers-Pehrson  used  a  simplified  2D  computer  code  to  model  the 
acceleration  of  the  metal  liner  by  the  explosive.  The  metal  was  assumed  to 
behave  as  a  fluid  and  its  thickness  was  ignored.  Each  point  along  the  metal 
surface  was  assumed  to  accelerate  according  to  the  equation 


V  = 


(1  -  exp  -  (t  -  t)/t) , 


(3.9) 


where  V  is  the  velocity  at  any  instant,  VQ  is  the  final  velocity,  T  is  the 
time  at  which  the  detonation  front  reaches  the  element,  and  x  is  the  time 
constant  for  the  acceleration.  VQ,  t  and  x  all  vary  with  initial  location. 
After  many  simulations  it  was  found  that  the  projection  angle  S  could  be 
described  by  the  expression 


V  V  'x 
►  o  o 

sinS  = 


(V  ’X) 
o 


(3.10) 


where  v  '  is  now  defined  as  being  the  derivative  of  V  with  respect  to 
distance  along  the  surface.  Typical  values  for  v  '  and  x  are 


.01  (mm/us)/mm  and  1  to  5  microseconds  respectively,  so  the  effect  of  these 
correction  terms  on  5  can  be  quite  significant. 


Chou  et  al  approached  the  problem  analytically  by  writing  the 
equations  of  motion  of  each  liner  segment  in  the  form 

(3.11 ) 
(3.1  2) 


6  =  ^  tan  (0  -  5), 

0  =  -  V*  cos  (0  -  6), 


where  0  is  the  angle  between  the  original  liner  contour  and  the  current  liner 
contour  at  1  (1  is  the  length  coordinate  along  the  liner  contour)  and  the  dot 
denotes  differentiation  with  respect  to  time.  Assuming  0-6  is  small, 
equations  (3.11)  and  (3.12)  yield 

6  =  -  I*  V’dt  +  ~  J*  (V2)'dt.  (3.13) 


When  equation  (3.9)  is  used  to  evaluate  this  expression  the  final  result 
becomes 


5  -  2§-It  V  +7T'  vo  (3-14> 

Note  that  the  first  correction  term  in  equation  (3.14)  is  the  same  as  found  by 
Randers-Pehrson.  Chou  et  al  have  compared  equation  (3.14),  (3.10)  and  (2.1) 
with  the  projection  angle  calculated  from  the  2D  Lagrangian  code  TEMPS.  The 
model  system  used  was  the  brl  81.3  mm  diameter,  42°,  conical  shaped  charge 
with  a  copper  liner  1 .9  mm  thick.  Fig.  7  shows  the  results.  The 
differences  between  the  two  non-steady  expressions  for  6  are  minimal,  and 
clearly  agree  much  better  with  the  TEMPS  calculation  than  the  original  Taylor 
expression. 


4.  CALCUIATION  OF  LINER  IMPUISE 


Before  the  non-steady  jet  formation  theory  described  in  the  previous 
section  can  be  used  to  model  e  real  shaped  charge  we  need  some  means  of 
calculating  the  velocity  imparted  to  the  liner  by  the  explosive  loading,  i.e. 
we  need  to  know  the  speed  VQ  and  projection  angle  6  as  a  function  of  length 
along  the  liner.  As  we  have  just  seen,  the  non-steady  Taylor  angle  formula 
provides  a  reasonably  accurate  relationship  between  these  two  variables,  so 
that  in  principle  we  need  determine  only  VQ  or  6.  Both  approaches  are  used 
in  shaped  charge  design. 


The  <?<5inmonly  used  method  of  tabulating  the  projection  angle 

for  a  given  charge  geometry  is  by  uSe  the  Richter  approximation  [23] 

as  modified  by  fc^oUbneaux  [24].  This  has  form 


1  1 


(4.1 ) 


where  K  and  <j>  aA  empirical  constants  which  ^apend  on  the  explosive  and  the 
angle  at  whicf>  tl<atonation  wave  intersects  the  liner.  They  can  be 
determined  frO^  experiments  with  explosive^ebl  slabs  or  sandwiches,  and  can 
also  be  found  hy  ^i^ting  to  a  small  subset  tf  the  shaped  charge  collapse  data, 
p  and  e  are  *HpSity  and  thickness  of  tb  liner  material  and  e  is  the 
explosive  thiC\n£^.  Harrison  [5]  has  modified  equation  (4.1)  to  take 
account  of  the  c<?hf£hement  of  the  casing  af^uM  the  charge.  He  writes  (4.1) 
in  the  form 


1_  1  KOe 

25  “  4>c  +  *e  ' 

where  B  is  a  C<^f^ti.on  factor  given  by 

8  »  1  +  A/tP^eo). 


(4.2) 


(4.3) 


p  and  e  are  Vs  density  and  thickness  of  tba  casing,  and  a  is  a  constant 
wRich  isCdeteA»ihe<i  from  experiment.  CarHoAe  and  Chou  [25]  have  also 
modified  equa^ph  (4.1)  to  take  into  account  ^he  non  planar  geometry  of 
conical  shaped  ^h^jr^es.  They  replaced  tbe  thickness  ratio  e/e  by  the  area 
ratio  A£/Ae,  ^  and  Ae  are  the  surface  a^as  generated  by  rotating  the 

lines  represeh^A^  hfje  thickness  e  and  e  ahiu't  the  cone  axis  for  each  point  x. 
An  example  of  accuracy  of  this  type  of  approach,  when  coupled  with  the  PER 
theory  of  jet  h?rmAion,  is  given  in  Fig.  8.  This  shows  a  comparison  of 
experimental  a^  theoretical  jet  velocity  distributions  for  the  heavily- 
confined  42°  shaped  charges  report^  hy  Di  Persio  et  al  [26].  It  is 

interesting  tj°  that  the  respectable  a^he^tneAt  between  theory  and 

experiment  is  tyAai-hed  using  only  the  steady  State  Taylor  equation  for  6,  and 
by  neglecting  Variation  of  K  and  <J>  wi^  ^he  angle  of  incidence  of  the 
detonation  wav^  c«n  hhe  liner.  The  values  were  K  =  0.25  cm^/g  and 

<j>  =23°. 


An  alt-^rh^tive  approach  to  the  pr^jAefl  of  calculating  the  liner 
impulse  is  to  use  a  form  of  Gurney  eq^Aiou  to  calculate  VQ  and  then  use 

one  of  the  aerations  "to  calculate  6*  The  original  Gurney  model  [18] 

(see  also  (27 ^  ^  discussion  of  this  work)  calculated  the  asymptotic  metal 

velocities  a^  a  ^uP^tion  of  the  explosive  t<j  nietel  mass  ratio  c/M  for 
symmetrical  <7enif^tAes  using  an  energy  equation  a.nd  the  assumption  of  a  linear 
velocity  grad^A  ih  the  explosive  gases,  a}0  could  also  treat  the  plane 
asymmetric  Including  a  momentum  ba^aA5?  equation.  The  formulae  have 

the  general 


V  //2E  =  f(C/M)  . 
o 


(4*4) 


where  /2E  is  an  empirical  constant  known  as  the  Gurney  energy,  and  the  form  of 
the  function  f  depends  on  the  geometry  of  the  explosive-metal  system.  This 
standard  Gurney  approach  is  unsuitable  for  modern  shaped  charge  design  in 
several  respects.  We  saw  in  the  previous  section  that  the  non-steady  Taylor 
angle  calculation  requires  that  the  full  velocity-time  curve  for  each  liner 
segment  be  known,  while  the  standard  Gurney  model  calculates  only  the 
asymptotic  velocity  VQ.  Also,  the  standard  model  applies  only  to  explosive 
systems,  and  needs  to  be  modified  to  handle  the  implosive  case.  Several 
papers  have  appeared  in  the  last  few  years  which  extend  the  Gurney  model  and 
attempt  to  remove  these  deficiencies. 

Chou  et  al  [16]  have  considered  the  problem  of  a  steady-state 
imploding  cylinder  using  the  same  general  assumptions  as  the  classical  Gurney 
theory.  Because  of  the  curved  geometry  in  the  imploding  case  the  momentum 
equation  contains  an  extra  term  which  represents  an  outward  impulse  applied  to 
the  explosive  gases.  This  extra  term  is  approximated  in  terms  of  the 
Chapman-Jouget  pressure,  the  characteristic  acceleration  time,  and  the  inner 
and  outer  radii  of  the  cylinder.  The  expression  for  the  final  liner  velocity 
as  a  function  of  charge-to-mass  ratio  is  then  compared  with  the  original 
Gurney  formula  and  a  2D  code  calculation,  and  is  found  to  be  a  significant 
improvement  on  the  standard  result. 

The  recent  paper  by  chanteret  [28]  extends  the  Gurney  method  in  the 
two  areas  needed  for  application  to  accurate  shaped  charge  calculations.  An 
analytical  model  is  developed  for  energy  partitioning  during  expansion  which 
allows  a  calculation  of  the  Gurney  energy  as  a  function  of  the  expansion 
ratio;  when  combined  with  the  Taylor  angle  equation  this  allows  the  entire 
velocity-time  curve  for  all  steady  state  geometries  to  be  calculated.  The 
results  are  found  to  agree  well  with  2D  code  calculations  or  with  available 
experimental  results.  chanteret  also  considers  the  problem  of  implosive 
geometries  by  introducing  the  concept  of  a  fictitious  rigid  boundary  in  the 
explosive  and  applying  this  to  the  case  of  the  imploding  cylinder  considered 
by  Chou  et  al.  A  comparison  of  liner  velocities  as  a  function  of  charge  to 
mass  ratio  c/M  shows  that  the  two  approaches  agree  remarkably  well. 

Hennequin  [29]  has  also  developed  an  analytical  model  of  liner 
collapse  based  on  assumptions  similar  to  those  of  Gurney.  The  laws  of 
conservation  of  mass,  momentum  and  energy  are  used  and  applied  in  detail  to 
the  explosive  gases,  liner  and  casing.  Both  cylindrical  and  conical 
geometries  are  considered. 

For  the  conical  case,  a  60°  apex  angle  was  used  and  liner  velocity 
as  a  function  of  time  for  different  elements  along  the  length  of  the  cone  was 
compared  with  2D  hydrocode  calculations  using  the  HEMP  code.  The  results 
agreed  well  for  the  first  1-2  ys,  but  then  showed  disagreements  of  the  order 
of  up  to  40%.  Nevertheless,  the  results  are  encouraging,  and  undoubtedly 
work  will  continue  in  this  direction. 

A  more  detailed  approach  to  the  calculation  of  liner  impulse  has 
been  made  by  Kiwan  &  Wisniewski  [30].  Their  idea  is  to  replace  the  liner  by 
a  system  of  discrete  solid  liner  elements  and  then  to  assume  that  both  the 
strength  properties  of  the  liner  material  and  the  interaction  forces  between 


liner  elements  can  be  neglected.  The  impulse  transmitted  by  the  detonation 
products  to  a  liner  element  is  then  obtained  by  integrating  the  excess 
pressure  along  the  path  of  the  liner  element.  The  pressure  wave  behind  the 
detonation  front  is  approximated  by  the  following  expression 


P(z)  =  P  -  P  +  (P  -  P  ) (z  -  1)/(1  -  X  ), 
f  o  f  o  o 

X  <  z  <  1  (4.5) 

o 

where  z  =  x/u_t,  P^  is  the  detonation  pressure,  PQ  the  ambient  atmospheric 
pressure  and  X  is  the  value  of  z  for  which  the  particle  velocity  vanishes  in 
the  Taylor,  wave.  Equation  (4.5)  applies  only  to  a  one  dimensional  plane 
Taylor  wave  in  a  tube  with  an  open  end  which  implies  a  completely  confined 
charge.  Kiwan  &  Wisniewski  show  how  this  can  be  modified  to  apply  to  an 
unconfined  charge,  and  also  to  take  account  of  the  cavity  in  the  explosive 
which  causes  a  reduction  in  the  total  available  energy.  The  liner  motion  is 
then  obtained  by  numerically  integrating  Newton's  equations  for  each  liner 
element,  i.e. 

a  xk  =  P(xk(t),  yk(t),  t)  sin  a  (4.6) 

a  yk  =  P(xk(t),  Yk(t) ,  t)  cos  a  (4.7) 


where  k  labels  the  liner  element,  P  is  given  by  an  equation  similar  to  (4.5), 
a  is  the  mass  per  unit  area  of  the  liner,  and  the  dot  indicates 
differentiation  with  respect  to  t.  The  jet  formation  equations  used  are 
similar  to  those  of  the  PER  theory  with  5  =0,  as  the  above  approach 
implicitly  assumes  that  each  element  of  the  liner  collapses  in  a  direction 
perpendicular  to  the  surface  of  its  initial  position  on  the  liner.  The 
theory  has  been  applied  to  two  shaped  charges,  having  diameters  of  3.44"  and 
3.60"  and  apex  angles  of  42°  and  60°  respectively.  Both  charges  are  loaded 
with  Comp.  B  and  have  .106"  copper  liners.  The  calculated  results  are  in 
good  agreement  with  the  experimental  results  of  Di  Persio,  Whiteford  and  Simon 
[31]. 

Kiwan  &  Wisniewski  [30]  have  remarked  that  a  possible  improvement 
for  their  model  would  be  to  compute  the  pressure  wave  using  one  of  the  many 
available  hydrocodes  and  then  to  couple  this  calculation  to  their  equations  of 
motion  for  the  liner  segments.  A  similar  approach  has  been  used  by  several 
groups  in  the  past  few  years  with  varying  degrees  of  success.  Edwards  & 
Godfrey  [32]  describe  a  method  based  on  the  CHAMP  computer  code  system  in 
which  the  high-explosive  detonation  and  hydrodynamics  are  computed  by  an 
explicit  Eulerian  finite  difference  method  while  the  metal  liner  is  viewed  as 
a  polygonal  line  of  mass  segments  in  which  the  metal  volume  and  internal 
stresses  are  neglected.  The  formation  of  the  jet  is  treated  analytically 
using  the  PER  theory.  This  approach  has  also  been  used  by  S.L.  and  H. 

Hancock  [33].  Van  Thiel  and  Levitan  [34]  describe  a  method  in  which  a 
Lagrange  code  (HEMP)  is  used  to  calculate  the  motion  of  the  detonation 
products  and  liner  segments,  which  are  then  coupled  to  analytic  expressions 


similar  to  those  of  PER  to  describe  the  jet  formation.  Kivity  et  al  [35] 
also  describe  a  scheme  in  which  the  explosive  and  liner  are  treated  by  the 
DISCO  code  (a  lagrange  code  similar  to  HEMP),  while  the  jet  formation  is 
described  by  the  PER  theory.  A  similar  approach  is  also  used  at  the 
Ballistic  Research  Laboratories  [36]  . 


5.  INVERSE  VELOCITY  GRADIENT  EFFECT  AND  JET  TIP  FORMATION 


The  PER  theory  is  based  on  the  assumption  of  a  negligible 
acceleration  time,  i.e.  it  is  assumed  that  the  detonation  impulse 
instantaneously  accelerates  the  liner  segment  to  its  terminal  velocity  v  . 

In  practice  this  does  not  occur.  We  saw  in  the  previous  section  that  each 
liner  segment  follows  a  velocity-time  curve  which  can  be  reasonably 
approximated  by  an  equation  of  the  form  of  (3.9),  where  the  time 
constant  t  can  have  a  value  of  several  microseconds.  This  means  that  liner 
segments  near  the  apex  of  the  cone  will  not  have  sufficient  time  to  be 
accelerated  to  their  terminal  velocity  VQ  before  colliding  on  the  cone  axis. 

As  the  distance  between  the  liner  segment  and  the  apex  of  the  cone  increases, 
more  space  (time)  is  available  for  the  segments  to  be  accelerated. 

Eventually  a  point  is  reached  at  which  the  segments  do  have  sufficient  time  to 
be  accelerated  to  v  ,  but  then  the  influence  of  the  decrease  of  the  explosive 
mass  to  liner  mass  takes  over  and  the  value  of  VQ  decreases  as  x  increases. 

The  point  at  which  this  usually  occurs  is  situated  about  one  third  of  the  way 
along  the  cone  axis  from  the  apex,  but  the  exact  position  depends  on  the 
geometry  used.  More  information  on  this  can  be  found  in  reference  [49] .  The 
net  result  is  that  the  velocities  of  the  liner  segments  just  before  collision 
with  the  axis  first  increase  to  a  maximum  value  somewhere  near  the  midpoint  of 
the  cone,  then  decrease.  The  jet  segments  produced  after  collision  of  the 
liner  segments  on  the  axis  will  also  follow  this  pattern.  A  further 
limitation  on  the  performance  of  the  shaped  charge  occurs  in  the  base  reqion 
where  the  combined  effect  of  a  significant  drop  in  C/M  ratio  and  the  explosive 
thickness  droping  below  the  critical  thickness  to  support  detonation  produce 
jet  elements  with  velocities  too  low  to  contribute  to  penetration  of  tarqet 
material.  This  is  described  more  fully  in  reference  [49] . 


The  velocity  distribution  within  the  jet  as  just  outlined  will  not 
be  stable  of  course.  In  the  leading  part  of  the  jet  a  particular  jet  element 
will  not  travel  for  long  before  the  jet  element  behind  it,  travelling  at  a 
higher  velocity,  will  collide  with  it  and  the  two  will  combine  to  form  a 
single  particle  moving  with  a  velocity  intermediate  between  the  velocities  of 
the  two  particles  before  collision.  This  process  of  collision  and  jet 
compression  will  lead  to  the  formation  of  a  relatively  massive  single  tip 
particle,  and  will  terminate  when  the  velocity  of  the  tip  particle  is  greater 
than  the  velocity  of  the  nearest  following  jet  segment.  This  phenomenon  has 
been  known  experimentally  for  some  time  [31],  and  the  presence  of  a  massive 
jet  tip  is  clearly  evident  in  flash  radiographs.  The  paper  by  ?erex  et  al 
provides  a  good  illustration  of  this  [37]  ,  as  does  the  MRL  flash  radioqraph 
reproduced  in  Fig.  2.  Note  that  the  high  quality  of  this  radiograph  also 
allows  the  ablation  of  the  jet  tip  particle  to  be  observed.  A  similar 
compression  occurs  for  the  slug,  as  has  been  noted  by  Kiwan  and  Wisniewski 
[30]  . 


Several  authors  have  successfully  modelled  the  formation  of  a 
compact  jet  tip  within  the  ID  approximation.  Carleone  et  al  [38]  have  made 
both  an  experimental  and  theoretical  study  of  this  effect.  They  describe  a 
series  of  experiments  using  the  81 .3  mm  BRL  Standard  Charge  in  which  the  apex 
end  of  the  liner  was  filled  to  different  heights  with  Woods  metal,  an  inert 
material  which  inhibits  the  collapse  of  the  apex  portion  of  the  liner. 

Several  radiographs  of  each  jet  were  taken  and  used  to  observe  the  jet  tip 
particle  and  to  measure  the  jet  velocity.  No  effort  was  made  to  directly 
examine  the  collapse  process  and  the  fate  of  the  woods  metal  however,  which 
could  effect  the  interpretation  of  the  data.  The  experimental  results  were 
compared  with  a  ID  theory  based  on  the  PER  equations,  but  containing  some 
important  modifications,  the  most  important  of  these  being  the  inclusion  of 
liner  acceleration.  The  liner  elements  were  still  assumed  to  travel  along  a 
straight  line  from  the  liner  to  the  axis,  but  each  element  was  assumed  to  have 
a  finite  constant  acceleration  until  it  either  hit  the  cone  axis  or  attained 
the  final  collapse  velocity  predicted  by  the  (steady  state)  Taylor  relation. 

It  was  assumed  that  the  jet  elements  experienced  a  perfectly  plastic  collision 
during  the  compression  stage,  and  the  velocity  of  the  tip  particle  was  then 
calculated  using  conservation  of  momentum.  Experimental  and  theoretical 
results  were  then  compared  and  showed  quite  good  agreement  for  jet  tip 
velocity  versus  percentage  of  ineffective  cone  height,  jet  velocity  versus  jet 
element  position,  and  jet  radius  versus  jet  segment  position. 


Harrison  [5]  has  described  an  iterative  scheme  for  the  calculation 
of  the  inverse  gradient  effect.  The  bending  angle  <{>  (which  is  twice  the 
Taylor  angle  6)  is  calculated  from  the  modified  Defourneaux  equation  (equation 
(4.2)).  A  new  equation  for  the  bending  angle  is  then  defined  as  follows: 


<j>  =  <j>e  (b/T*) 

N 


V2 


(5.1 ) 


where  T*  is  the  time  which  an  element  takes  to  reach  the  axis,  e  is  the 
thickness  of  the  explosive,  and  b  is  related  to  the  density  and  thickness  of 
both  the  liner  and  explosive.  The  iteration  is  carried  out  between  equation 
(5.1)  and  the  expression  for  the  collapse  time  T*  until  <{>  approaches  <{>  to 
within  some  desired  degree  of  accuracy.  The  scheme  thenNproceeds  in  a 
similar  manner  to  the  calculation  described  by  Carleone  et  al  [38] .  Harrison 
makes  an  extensive  comparison  between  this  simple  ID  theory,  a  2D  hydrocode 
(the  BRL  version  of  the  HELP  code)  and  experiment  for  a  variety  of  different 
shaped  charges  and  generally  finds  quite  good  agreement. 


Kiwan  and  Wisniewski  [30]  have  also  made  a  detailed  calculation  of 
the  "pile  up"  of  the  leading  jet  elements  to  form  a  jet  tip  particle,  but 
their  comparison  with  other  calculations,  or  with  experiment,  is  nowhere  near 
as  extensive  as  either  Carleone  et  al  or  Harrison.  One  novel  feature  of 
their  calculation  is  that  they  allow  the  collisions  in  the  jet  compression 
region  to  be  either  purely  elastic  or  purely  plastic.  This  makes  a 
difference  of  approximately  10%  co  the  final  jet  tip  velocity. 
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6.  EFFECT  OF  COMPRESSIBILITY  ON  JET  FORMATION 


The  original  steady  state  shaped  charge  theory  of  Birkhoff  et  al, 
and  its  extension  to  the  non  steady  case  by  Pugh  et  al,  both  assumed  that  the 
liner  material  flowed  as  an  incompressible  fluid.  The  result  of  this 
assumption  of  incompressibility  is  that  a  coherent  jet  must  always  be 
formed.  If  the  fluid  is  compressible  however  the  situation  is  more 
complicated,  and  there  is  the  possibility  of  either  a  no  jet  configuration,  or 
the  formation  of  a  non  cohesive  jet  [39]. 


The  2D  wedge  shaped  geometry  formed  by  the  angular  impact  of  two 

plates  at  high  velocity  has  been  analysed  in  detail  by  Walsh,  Shreffler  and 

Willig  [39],  and  by  Cowan  and  Holtzman  [40].  They  find  that  if  the  fluid 

velocity  of  the  collision  is  subsonic  (relative  to  a  frame  moving  with  the 

velocity  of  the  collision  point)  then  a  coherent  jet  will  always  be  formed,  as 

in  the  incompressible  case.  For  supersonic  collisions  however  jetting  only 

occurs  if  0  >  0  ,  where  0  is  the  critical  angle  at  which  the  shock  wave  at 
c  c 

the  stagnation  point  becomes  detached  (the  angle  between  the  plates  is  20). 

For  0  <  0  no  jetting  occurs, 
c 


These  conclusions  have  been  verified  by  Harlow  and  pracht  [41 ]  by 
numerical  calculations  of  the  oblique  collapse  of  two  metal  plates  using  the 
Particle-in-Cell  method.  These  calculations  were  also  able  to  include  the 
effect  of  viscosity  on  the  collapse  process  and  were  also  able  to  examine  the 
transient  stages  leading  to  the  formation  of  the  steady  state. 


Chou,  Car leone  and  Karpp  [42,43]  have  reviewed  this  earlier  work  and 
considered  its  application  to  the  axially  symmetric  geometries  which  occur  in 
conically  lined  shaped  charges.  They  note  that  the  criterion  given  by  Walsh 
et  al  and  by  Cowan  and  Holtzman  is  based  on  the  maximum  angle  for  an  attached 
shock  at  the  collision  point  as  determined  by  the  shock  polar  for  the 
material.  As  the  shock  conditions  and  shock  polar  are  applicable  to  curved 
shocks  at  individual  points  on  the  shock,  they  therefore  believe  that  the 
above  jetting  criterion  should  also  be  applicable  to  a  shaped  charge  with  a 
conical  liner.  They  therefore  propose  the  following: 


Jetting  Hypothesis 

1 .  For  subsonic  collisions,  a  solid  coherent  jet  always  occurs. 

2.  For  supersonic  collisions,  jetting  occurs  if  0  >  0  for  a  given  v; 
the  jet  formed  is  not  coherent. 

3.  For  supersonic  collisions,  if  0  <  0  for  a  given  v,  no  jet  will  be 

c 

formed. 

(Here  v  is  the  velocity  of  the  liner  in  the  frame  moving  with  the  velocity  of 
the  collapse  point).  The  remainder  of  the  report  by  Chou  et  al  [42]  details 
a  large  number  of  experiments  and  numerical  computations  which  are  in 
agreement  with  their  hypothesis. 
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A  similar  jetting  hypothesis  has  been  used  by  Harrison  in  the  BASC 
code  [5],  His  criterion  for  a  coherent  jet  is  that  the  velocity  v  should 
satisfy  the  inequality 


v  <  1  .23  C  (6.1  ) 

where  C  is  the  sound  velocity  in  the  liner  material  and  the  numerical  factor 
has  been  determined  by  comparing  jet  radiographs  with  BASC  results. 


Finally,  we  note  that  an  interesting  application  of  this  work  has 
recently  been  made  by  Carleone  and  Chou  [44]  in  calculating  a  theoretical 
maximum  velocity  for  a  coherent  shaped  charge  jet.  Using  the  above  jetting 
hypothesis,  the  simple  Gurney  model,  and  some  geometrical  arguments,  they 
derive  the  following  upper  limit  for  jet  velocity 

V.  <  2.41  C  (6.2) 

3 

Carleone  and  Chou  compare  this  criterion  with  a  large  amount  of  published  and 
unpublished  shaped  charge  data  and  conclude  that  the  maximum  jet  tip  speed 
given  by  equation  (6.2)  has  not  been  exceeded  by  any  cohesive  jet. 

It  is  interesting  to  note  that  the  application  of  equation  (6.2)  to 
a  steel  liner  results  in  a  theoretical  V-  ax  of  9.8  mm/ps  (the  value  for 
copper  is  10.2mm/jis),  whereas  Mader  et  al  have  recently  reported  the  formation 
of  a  jet  from  a  4.0  mm  thick  hemishell  having  a  jet  trip  velocity  of 
18mmus  [48].  it  is  uncertain  however  whether  the  theory  of  Carleone  and  Chou 
is  applicable  to  this  geometry. 
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GLOSSARY  OF  SYMBOLS 


Velocity  of  detonation  wave 

Final  collapse  velocity  of  liner 

Original  cone  half  angle 

Generic  symbol  for  collapse  angle 

Collapse  angle  in  Steady  State  theory 

U  /cos  c 
D 

Taylor  angle 

Velocity  of  stagnation  point 

Velocity  of  liner  in  frame  moving  with  velocity  of  stagnation  point 

Speed  of  jet 
Speed  of  slug 

Mass  per  unit  length  of  liner 

Mass  per  unit  length  of  jet  at  stagnation  point 
Mass  per  unit  length  of  slug  at  stagnation  point 
Measures  distance  along  cone  axis  from  apex 
Cylindrical  coordinates  of  typical  liner  element 
Time  for  detonation  wave  to  travel  distance  x 
a  +  f 

Derivative  of  VQ  with  respect  to  x  or  l ,  where  appropriate 
Difference  between  0  and  0+ 

Liner  velocity  when  finite  acceleration  time  is  considered 
Time  constant  for  acceleration  of  liner 

Angle  between  original  liner  contour  and  current  liner  contour 
Length  corrodinate  along  liner  contour 
Derivative  of  8  with  respect  to  time 


Empirical  constants  used  in  Richter  approximation  for  5 
Density  of  liner  material 
Thickness  of  liner  material 
Thickness  of  explosive 

Correction  factory  for  Richter  approximation  to  take  casing  into 
account 

Density  of  casing  material 
Thickness  of  casing  material 

Experimentalloy  determined  constant  for  use  in  modified  Richter 
approximation 

Surface  areas  generated  by  rotating  the  lines  representing 
thicknesses  e  and  e  about  the  cone  axis 

Explosive  to  metal  mass  ratio 

Gurney  energy 

Pressure  wave  behind  detonation  front 
Deontation  pressure 
Atmospheric  pressure 

Scaled  coordinate  along  cone  axis,  equal  to  x/UDt 

Value  of  z  for  which  particle  velocity  vanishes  in  Taylor  wave 

Mass  per  unit  area  of  liner  material 

Cartesian  coordinates  of  typical  liner  element 

height  of  cone 

Bending  angle  (twice  Taylor  angle) 

Nth  iteration  of  bending  angle 
Constant  used  in  bending  angle  formula 
Time  taken  for  liner  element  to  hit  axis 

Critical  angle  at  which  the  shock  wave  at  the  stagnation  point 
becomes  detached 


Velocity  of  sound  in  liner  material 


* 
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Figure  4.  Diagram  for  determining  Taylor  angle  formula  (from  reference  12) 


